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1.0 Introduction 



For decades the U.S. Navy has used large steel barges as a vital part of their 
Military Preposition Force (MPF) offload process. The MPF is a fleet of cargo ships 
prepositioned at strategic locations ready for deployment when called upon. For any 
extended military action, the MPF is the main source of equipment transport. Many 
of the locations to which the U S military may be called into action do not have deep 
water facilities that can accommodate these vessels. As a result, the equipment is 
offloaded several miles off shore and transported to the beach or permanent/temporary 
pier via barges or "causeways" as the Navy has designated them. 

The causeways currently being used by the Navy measure 21 foot (ft) wide, 90 
ft long and have a 5 ft draft. The shape is very close to a rectangular steel box. Due 
to it’s shallow draft and particular geometry, the causeways quite often have water on 
the deck even in relatively low sea states. Hence their stability is sensitive to wave 
action. To improve the reliability of their operation, a new generation of causeways, 
in which the dimensions increase to 25 ft wide, 120 ft long and a draft of 8 ft is 
currently under conceptual design and is the focus of this analysis. 

1.1 Objective 

A series of experiments have been conducted to provide data to calibrate the 
prediction capabilities of several time, frequency and probability domain computer 
programs developed by the Naval Facilities Engineering Service Center (NFESC), 
Amphibious Division. Of particular interest to the NFESC engineers is a study of 
capsizing when the causeways are subjected to high sea states for relatively long 
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durations. Simple low degree-of-freedom dynamic models of the barge motions have 
been developed for stochastic analytical and numerical studies to estimate the 
probability of capsize. The motivation behind this study is to identify suitable system 
damping parameters for the single-degree-of-freedom (SDOF) model. 

1.2 Literature Review 

The motion of free-floating vessels in open seas often contain complicated 
nonlinear behaviors due to large displacement and hydrodynamic effects. In 
particular, the stability of roll motion of barges with high centers of gravity is of 
great practical importance. Existing stability criteria are expressed in terms of 
minimum values of certain key features of the righting arm. For certain classes of 
vessels, static stability standards based on statistical and other analyses of intact static 
condition are sufficient for design purpose, and can give a qualitative understanding of 
the stability behavior for the naval architect (Soliman and Thompson 1991). Although 
this curve is found to be an important vessel characteristic in assuring safety 
(Falzarano et al. 1992), other vessel properties are also significant. These include 
hydrodynamic and viscous roll damping, as affected by the wave exciting force. 

Because the cargo is carried above deck on a causeway, the center of gravity 
of the entire system is considerably higher in proportion to that of a conventional ship 
which carries most of its cargo below deck. The dynamics of a high center-of-gravity 
barge rolling in seaways can become difficult to model and predict because of the 
highly nonlinear characteristics encountered. Accidental capsize of log barges (which 
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have very similar geometry and loading conditions as those considered in this study) 
under average to calm sea condition has been well documented (McAllister 1995). 

1.3 Scope 

The scope of this study is to model the barge response when subjected to beam 
sea conditions as a SDOF system for roll motion. The measured wave, either regular 
or random, is input into a computer program which uses a fourth-order Runge-Kutta 
integration method to numerically predict roll motion response. The simulated results 
are compared to measured data to determine the best system damping parameters. 

Four different forms of the damping moment of the SDOF model are analyzed. In 
each form all coefficients are known except for the linear and nonlinear damping 
parameters. Each form uses a combination of relative motion "Morison" damping 
and/or linear structural damping. Predicted results of the four forms of the damping 
moment are examined to determine the most suitable model. A sensitivity study on 
the response to various system parameters is then conducted on the selected form. 

Each form of the model uses a thirteenth-order polynomial restoring moment. An 
analysis is also conducted to determine the effects of using lower order terms to 
represent the restoring moment. Finally, results of the SDOF model from this report 
are compared with a 3DOF and a 2DOF model subjected to the same wave 
environment to determine which model more accurately predicts roll response. 



4 



2.0 System Description 

The system considered in this analysis is that of a free floating barge subjected 
to regular and random waves in open seaways. In general, a six degree-of-freedom 
(6DOF) model is required to fully specify the motions of the barge. Subjecting a 
barge to beam seas would be considered the most critical scenario for barge capsize. 
This simplifies the equations of motion because under ideal beam conditions, the 
response motions for a barge may be described by a three degree-of-freedom (3DOF) 
model with heave, sway and roll represented by 

my + m^cos(|^)(y-v) + m a j in(||^|)(y-i>) + C n (y-v) 

+ C 22 K 0 i-v )l^" v l ” ~ m(z g cos<t>)4> 

+ R 33 (<j>,z,v,^L)sm(fy + Kmoor y = 0 

mz + m a cos + m sin(|^|)(z->v) + C 33 (z-w) 

" dy ^ dy 1 (2.1) 

+ C 33 (z-w)| z-w| + m4>y - m(z g cos<j>) 4> 2 + mg 

+ R 33 (z,<i>, cos A = 0 

dy dy 

lj * i.ji- 15) . c (i- ii) . c I 

dy L dy s dy dy 
+ m (z g cos<t>)(t>z - m(z g cos<f>)y + z ,t],^)cos(^) 

- mgz g sin<f> = 0 

(see Bartel 1996 for derivation and definitions). Although this model is significantly 
simpler than a large number of time domain large body analysis models based on 
potential theory, it is evident that these equations are still relatively involved and 
computationally intensive to solve numerically compared to that of a SDOF model. 
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The 3D0F model based on these equations is currently being analyzed under a 
separate study (Bartel 1996). 

2.1 Assumptions 

There are several assumptions made in the development and application of the 
model equations for this analysis. The first is that the nominal wave length is 
significantly greater than the beam of the barge. As a result, the wave surface can be 
assumed linear across the beam. Secondly, wave forces and moments act at the 
center of gravity in lieu of integrating pressure over the surface. Finally, it is 
assumed that roll response can be uncoupled from heave and sway, therefore the 
focus will be only on the pure roll motion of the barge modeled by a SDOF system. 

2.2 Equation of Motion 

Uncoupled roll motion subjected to beam sea conditions, assuming long wave 
excitation and deep water conditions, is modeled by the equation 

lj> * I^ti-dy/dy) + C^-dijIdy) + C MN (4>-dr } /dy)\j>-dy/dy\ (2.2) 

+ RJ,4>,ri,drildy) = 0 

where I M is the rotational inertia (in air) of the vessel about the roll axis, I A44 the 
hydrodynamic added inertia, <j> the roll acceleration, dq/dy the wave slope 
acceleration, C UL the linear damping coefficient, 4> the roll velocity, dq/dy the wave 
slope velocity, the nonlinear damping coefficient, and R u (4>,q,dq/dy) the 
nonlinear restoring moment. 
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2.3 Description of Wave Field 

Based on linear wave theory, the velocity potential and related wave profile, 
time derivatives and partial spatial derivatives are as follows: 



* = A± eo **Q , **> wy - «o 

co cosh (Jdi) 

V = = Acos(ky - cor) 

g dt 

n = — = co/1 sin (ky - cor) 
dt 

ff = = -co 2 /lcos(;ty-cor) 

at 

i)' = — =- £/lsin(fcy - cor) 

5y 

ij' = = o)kAcos(ky - cor) 

3y 

i)' = — = co 2 £v4sin(£y - cor) 

ay 

If deep water wave conditions are applied and water particle characteristics are 
considered at the surface, the wave number, k, can be represented by 



(2.3) 



k= 



CO 



g 



(2.4) 
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and the resulting wave field characteristics reduce to 

ij = /I cos (Ary - <j)t) 

77 = avlsin (Icy - a it) 

7 

= —co 77 

= -f * 0-5) 

o > 3 

= —V 

g 

a> 3 . 

= —V 

g 

2.4 Inertia and Added Inertia Coefficients 

For a given barge, the rotational inertia, 1 H , can be analytically computed. 

The hydrodynamic added inertia, 7^, depends on the frequency of the barge roll 
motion. The value of the added inertia under the range of frequencies of the roll 
motion response observed in this study is relatively constant. Therefore, a constant 
added inertia coefficient is adopted for convenience of analysis. The first two terms 
of equation 2.2 represent the inertial moment of the barge in air, I u 4>, and the relative 
motion added inertia moment, I AU (4>-dfi/dy), respectively. The term ( <j>-dfi/dy ) is the 
relative motion between the roll acceleration of the barge and the wave slope 
acceleration resulting from wave excitation. For this analysis, the value of 1^ = 
2.161046e6 slug-ft 2 (Paulling 1995). With this information and knowing the natural 
period of the barge is approximately 5 seconds, the added inertia is calculated to be 
1^44 = 1.3725e6 slug-ft 2 . 



V 
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2.5 Damping Moment 

In equation 2.2, the damping moment resulting from radiation of waves from 
the barge (corresponding to linear energy dissipation in potential theory) is 
represented by the linear term, C Ml £4>-di]/dy). The effects of damping resulting from 
flow separation and turbulence are represented by the nonlinear Morison term, 

Q 4 nW> • drj/dy) | 4> - dij/dy \ . C UL and are the linear and nonlinear damping 
coefficients expressed in the following convenient form: 

r = 2 * w „*£l ^ = 2*u n *| ;v (2.6) 

C 44i~ r U 44V J 

l T l T 

where co„ is the "natural frequency" of the barge (determined by small amplitude free 
vibration tests) and is equal to approximately 5 seconds (Yim et al. 1995). % N and | L 
are the damping parameters associated with the linear and nonlinear damping terms 
(to be determined by system identification and/or numerical iteration). The terms 
have been normalized by the sum of the inertia and added inertia terms, I T =l u + 

T^. The term (4> - dij/dy ) is the relative motion between the roll velocity of the barge 
and the wave slope velocity resulting from wave excitation. 

2.6 Roll Restoring Moment 

The restoring moment term, R^^^^rj/dy), represents the moment created 
from buoyancy and gravity forces acting on the barge under static conditions. The 
restoring moment is a thirteenth-order (13th-order) polynomial made up of 
coefficients of the righting moment curve multiplied by the corresponding power of 



the relative motion between the roll displacement and the wave slope as shown in the 
following equation: 
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p /JL „ a /a „, B l3 ' 5 (<i>-dr,/dy)\ B u 5 (4> -drj/dy) 3 B 9 ,(<t>-d V /dy) s 

B 44 (<t>,y,ai]/dy) = j + - + - + 

£ 75 (<MVd)0 7 B^tf-dv/dy) 9 B 3 s ( 4> -dy/dy) 11 B ls (4>-dr,ldy) a ^ 

’ - - - ■ + + I + - 

if 1 

The coefficients B, 5 , B 3 5 , etc. are established by fitting a thirteenth-order polynomial 
to the righting-moment curve shown in Figure 2.1 (Yim et al. 1995). The term is 
normalized by the total inertia, I T . Although a lower-order polynomial can be used 
with sufficient accuracy for small angle motions, the additional time and effort 
associated with using a thirteenth-order polynomial is minimal, and in the case of 
large angle motions, the higher order terms are essential. 

2.7 Various Forms of Roll Equation of Motion 

Several combinations of the terms representing the damping moment in the 
equation of motion are examined here. In addition to equation 2.2 which uses relative 
motions in both the linear and nonlinear damping terms, models which approximately 
account for and neglect different physical effects will be analyzed. One combination 
assumes the magnitude of the wave-slope velocity is negligible compared to the roll 
velocity (e.g. near roll resonance). As a result only the roll velocity component will 
be considered in the linear and nonlinear damping terms. This is modeled by 
equation 2 . 8 : 
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Figure 2.1 13th Order Roll Restoring Moment Curve. - numerical (Paulling 1995) 

-- polynomial fit (Yim et al. 1995). B 135 = 5.6051e6, B u5 =-1.1657e6, 

B,^ 1.7258e6, B^ 9.6915e6, B 5i5 = -9.3298e6, B 3i5 = 3.3999e6, B 1i5 = -4.4569e6. 
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lj> + I^ti-dn/dy) + + C^l^l (2.8) 

+ R^^^ldy) = 0 

Equation 2.9 below assumes that the linear damping coefficient is dominated 
by the barge characteristics (damping induced by the presence of mooring lines in the 
experiments which is not modeled explicitly in this study). The relative motion 
between the barge and wave slope velocity is accounted for only in the nonlinear 
damping term. 



J m 4> + ^ 444 -dfj/dy) + C ml 4> + C MN (4> -dy/dy) | 4> -dij/dy \ 

+ R A4 (<t>,'n,d , n/dy) = 0 

The final form of the model to be analyzed is equation 2. 10: 



(2.9) 



l u 4> + l^ty-dri/dy) + C mn ( 4> -drj/dy) | 4>-di]/dy \ 
+ Rm = 0 



( 2 . 10 ) 



This model assumes that damping effects are strictly nonlinear and relative motion 
between the wave slope velocity and barge roll velocity significantly influences the 
roll response. 
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3.0 Model Tests and Identification of Damping Parameters 

The four different forms of the model for roll motions are analyzed to 
determine which form best predicts roll response when the barge is subjected to beam 
sea conditions. In each case, the prototype measured wave is sampled at 2 Hertz 
(Hz). The data is then filtered so that frequencies above 0.25 Hz are disregarded 
because they do not contribute to the response motion. Wave profile and water 
particle velocity as well as wave slope, wave slope velocity and wave slope 
acceleration are then numerically computed. Regular wave cases analyzed are that of 
the 5, 6, 8, and 10 sec waves with a desired 6 ft wave height. In addition, two 
random cases are modeled; a Bretschneider spectrum of 6 ft significant wave height 
and 8 sec dominant wave period, and a white noise spectrum of 6 ft significant wave 
height and wave periods ranging from 4 to 20 secs. 

In each of the forms, all terms are known except for the linear and nonlinear 
damping parameters, and % N used in the damping coefficients, C ML and C^. An 
initial estimate for the damping parameters is obtained through system identification 
procedures. By comparing the standard deviation of predicted roll response to 
measured roll response, an iterative approach is then used to identify the most 
appropriate damping parameters for the model. The inertia coefficients and the 
nonlinear restoring moment coefficients, constant for all the wave cases, are 
determined through experimental measurements or computer analysis. The other 
coefficients are determined through the measured data for each particular wave 



excitation test. 
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3.1 Model Test Procedure 

Regular and random wave model tests have been conducted on a single barge 
to simulate response in open seas. The facility used to conduct the tests was a 48 ft 
wide, 250 ft long basin which was capable of generating long crested waves up to 3 ft 
in height. A 1/16 scale model of the 120 ft long, 25 ft wide and 8 ft deep full scale 
barge was used. For this analysis, only results of the barge subjected to beam sea 
conditions are examined. Although the objective of the tests is to simulate motions of 
the barge in open seas, mooring lines were attached to the model to prevent it from 
drifting down the basin. The length, weight and sag of the mooring lines are 
designed to ensure the natural period of the mooring system would be much larger 
than the wave excitation period. 

The particular wave conditions generated for this report were 4 regular wave 
and 2 random wave tests. Regular waves of approximately 6 ft wave height (full 
scale) were conducted for wave periods of 5, 6, 8 and 10 secs. The random wave 
tests were for a Bretschneider spectrum of 6 ft significant wave height with and 8 sec 
dominant wave period, and for a band limited "white noise" test of 6 ft significant 
wave height with wave periods ranging from 4 to 20 secs. Model tests were run until 
steady state had been reached before data collection began. Only 1024 secs (17 
minutes) of data was collected for each run in order to avoid significant energy 
accumulation from wave reflection, re-reflection and other "tank noise" which might 
build up over time and could corrupt the wave excitation. 
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3.2 Description of Analytical Prediction Procedure 

Test SB26, the 6 ft, 5 sec regular wave, will be used as an example of the 
procedure followed to establish damping coefficients and simulated roll response. 

The entire time series of the measured wave input is examined to confirm the desired 
wave height and period was consistent with the target spectrum. Figure 3.1a shows 
the 1024 sec time series used in this analysis. The measured data is then subjected to 
a low pass filter with cutoff frequency of 0.25 Hz after digital sampling. The filtered 
wave is then displayed as shown in Figure 3.1b, this is the wave data input into the 
simulation program. To ensure the input wave for simulation is correct, the spectral 
densities of the measured and filtered wave are compared as shown in Figure 3.2. 
The damping parameters for the linear and nonlinear damping coefficients are input 
into the model equation and the simulation is executed. The results of the measured 
roll response, Figure 3.1c, are then compared to those of the simulated roll response, 
Figure 3. Id. A small section of the time series is compared as well, as shown in 
Figure 3.3. The standard deviation of the measured and simulated data is then 
compared to ensure the most accurate match of roll response. In the random wave 
tests, the above procedure is performed and an additional comparison of the spectral 
and probability densities of roll response are made to ensure the most accurate match 
of roll response. 



roll (degrees) roll (degrees) 
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(a) Measured Wave Profile 




(b) Filtered Wave Profile 
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(d) Predicted Roll time (seconds) 

Figure 3.1 Comparison of Measured and Filtered Wave, Measured and Predicted 
Roll, Regular Wave H= 6.6 Ft, T= 5 Sec (SB26) 
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Figure 3.2 Comparison of Spectral Densities (Ft^2/Hz), Measured and Filtered Wave 
(SB26), - measured, — predicted 
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time (seconds) 



Figure 3.3 Comparison of Measured and Predicted Wave and Roll, Regular Wave, 
H= 6.6 Ft, T= 5 Secs (SB26) 

- measured wave, predicted wave, - measured roll, .. predicted roll. 
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3.3 Form 1: Relative Motion Damping Moment 

The first model to be analyzed is that of equation 2.2: 



J^4> + Ijuit-drj/dy) + C UL (<t>-dri/dy) + C MN (j> -dy/dy) \ <i> -drj/dy \ n 2) 

+ Ru(4>,V,tydy) = 0 



Recall this form assumes both the linear (radiation damping) and nonlinear (Morison 
type local flow separation) have significant influence on the roll motion response. 
The resulting linear and nonlinear damping coefficients for this model are listed in 
Table 3.1: 



Wave Test 


Linear Damping 
Parameter 


Nonlinear Damping 
Parameter 


H=6.6 ft, T=5 sec 
, (SB26) 


0.19 


0.19 


H=6 ft, T =6 sec 
(SB27) 


0.11 


0.11 


H=7.2 ft, T=8 sec 
(SB29) 


0.80 


0.20 


H=5.7 ft, T=10 sec 
(SB30) 


0.03 


0.07 


Bretschneider 
Spectrum (SB25) 


0.04 


0.04 


White noise 
Spectrum (SB33) 


0.08 


0.06 



Table 3.1 Linear and Nonlinear Damping Parameters for Form 1, Equation 2.2. 
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3.3.1 Regular Wave, H= 6 Ft, T= 5 Sec (SB26) 

Figure 3.4 represents the roll response of the model barge subjected to a 6.6 
ft wave with a 5 sec period. The maximum roll response is approximately 13 
degrees. The model indicates a linear damping parameter of 0.19 and a nonlinear 
damping parameter of 0.19. The amplitude of the measured and predicted response is 
nearly identical. However, the predicted roll response is out of phase by 
approximately 1 sec (1/5 period) with the measured data. 

3.3.2 Regular Wave, H= 6 Ft, T= 6 Sec (SB27) 

Figure 3.5 represents the results of the model barge subjected to a 6 ft regular 
wave with a 6 sec wave period. The maximum roll response is 11.5 degrees. The 
model indicates both the linear and nonlinear damping parameters are 0. 11. Figure 
3.5 shows that amplitude of the predicted roll response matches the measured data 
well. As in the 5 sec wave case, the predicted roll response is approximately 1 sec 
(1/6 period) out of phase with the measured data. 

3.3.3 Regular Wave, H= 7.2 Ft, T= 8 Sec (SB29) 

Figure 3.6 represents the results of the model barge subjected to a 7.2 ft 
regular wave with an 8 sec wave period. The damping parameters established for this 
wave test are 0.80 for linear damping and 0.20 for nonlinear damping. Note that the 
linear damping parameter is significantly higher than the two previous test cases. The 
maximum roll response of the measured data is 4. 1 degrees and is well matched by 



roll (degrees) and wave (feet) 
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Figure 3.4 Comparison of Measured and Filtered Wave, Measured and Predicted 
Roll, Regular Wave, H= 6.6 Ft, T= 5 Sec (SB26, eqn 2.2) 

- measured wave, -. predicted wave, — measured roll, .. predicted roll. 
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time (seconds) 



Figure 3.5 Comparison of Measured and Filtered Wave, Measured and Predicted 
Roll, Regular Wave, H= 6 Ft, T= 6 Sec (SB27, eqn 2.2) 

- measured wave, -. predicted wave, — measured roll, .. predicted roll. 
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time (seconds) 



Figure 3.6 Comparison of Measured and Filtered Wave, Measured and Predicted 
RoH, Regular Wave H= 7.2 Ft, T= 8 Sec (SB29, eqn 2.2) 

- measured wave, -. predicted wave, — measured roll, .. predicted roll. 
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the predicted response of the model. This form of the model is not able to 
consistently match the phase of the response. The predicted data appears to be 
approximately 1 sec (1/8 period) out of phase (at the peak) with the measured data. 

3.3.4 Regular Wave, H= 5.7 Ft, T= 10 Sec (SB30) 

Test SB30, in which the excitation is a 5.7 ft wave with a 10 sec period, proves 
to be the most difficult response to match by any form of the model. Because the 
natural period of the barge is approximately 5 secs, which is a multiple of the wave 
excitation period, the response contains a superharmonic. 

Figure 3.7 represents the best comparison of standard deviation and 
superharmonic response which could be obtained with the given form of the model 
equation (equation 2.2). If the damping parameter is increased, the superharmonic 
characteristics are eliminated. If the damping parameter is decreased, the simulated 
response amplitude is much greater than the measured response amplitude. The 
damping parameters used in this figure are 0.03 for the linear coefficient and 0.07 for 
the nonlinear coefficient. The standard deviation for the measured response is 
calculated to be 1.84. The standard deviation of the simulated response is 1.90. 

Figure 3.8 shows a good comparison of the measured and simulated spectral densities 
for roll response. 



wave (feet) and roll (degrees) 
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Figure 3.7 Comparison of Measured and Predicted Wave, Measured and Predicted 
Roll, Regular Wave, H= 5.7 Ft, T= 10 Sec (SB30, eqn 2.2) 

- measured wave, -. predicted wave, — measured roll, .. predicted roll. 
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Figure 3.8 Comparison of Spectral Densities (Deg^2/Hz), Measured and Predicted 
Roll, H= 5.7 Ft, T= 10 Sec (SB30, eqn 2.2), - measured, — predicted 




Figure 3.9 1024 Sec Time Series of Predicted Roll, H= 5.7 Ft, T= 10 Sec 
(SB30, eqn 2.2) 
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Figure 3.9 provides a representation of predicted roll response for the entire 1024 sec 
time series. It is evident that the superharmonic does not occur in all locations with 
the same amplitude. A parametric study presented in a later section of this report will 
provide a better understanding of the effects of varying the linear and nonlinear 
damping parameters. 

3.3.5 Bretschneider Spectrum, Hs= 6 Ft, Tp= 8 Sec (SB25) 

Figure 3.10a - 3.10d represent the measured and filtered wave input and the 
resulting measured roll response and predicted roll response (to filtered wave 
excitation) for a Bretschneider wave spectrum. As shown in Figure 3.11, the spectral 
densities of the measured and filtered wave excitation input are compared to ensure 
the proper input is being used for prediction. The damping parameters used for 
prediction are equal to 0.04 for both the linear and the nonlinear terms. A visual 
comparison of Figures 3.10c and 3.10d indicates the model provides an accurate 
prediction. The standard deviation for the predicted roll response is 3.95 which 
compared favorably to the standard deviation of 3.95 for the measured data. Figure 
3.12 indicates the spectral densities are a very close match. The final comparison 
made between measured and simulated data is a histogram of roll response. Figures 
3.13a and b confirm that the model provides an accurate representation of roll 
response motion. 



roll (degrees) roll (degrees) 



27 



5 




5 1 1 1 1 1 i_i 

0 200 400 600 800 1000 



(a) Measured Wave Profile 
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(b) Filtered Wave Profile 
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(c) Measured Roll 




0 200 400 600 800 1000 

(d) Predicted Roll time (seconds) 

Figure 3. 10 Comparison of Measured and Filtered Wave, Measured and Predicted 
Roll, Bretschneider Spectrum, Hs= 6 Ft, Tp= 8 Sec (SB25, eqn 2.2) 
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Figure 3.11 Comparison of Spectral Densities (Ft A 2/Hz), Measured and Filtered 
Wave, Bretschneider Spectrum Hs= 6 Ft, Tp= 8 Sec (SB25, eqn 2.2) 

- measured, — predicted 
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Figure 3. 12 Comparison of Spectral Densities, Measured and Predicted Roll 
(Deg^2/Hz), Bretschneider Spectrum , Hs= 6 Ft, Tp= 8 Sec (SB25, eqn 2.2) 
- measured, — predicted 



occurance occurance 
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(a) Measured Roll Amplitude (degrees) 




Figure 3.13 Histograms of Measured and Predicted Roll, Bretschneider Spectrum 
(SB25, eqn 2.2) 
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3.3.6 White Noise Excitation, Hs= 6 Ft (SB33) 

Figures 3.14a through 3.14d represent the measured and predicted wave input 
and measured and predicted roll response for a white noise excitation with wave 
periods of 4 to 20 secs and a 6 ft significant wave height. A comparison of the 
spectral densities, Figure 3.15, confirms that the simulated wave input matches the 
measured wave excitation. The linear and nonlinear damping parameters used in the 
model are 0.08 and 0.06 respectively. A visual comparison of Figures 3.14c and 
3. 14d indicate the model provides a good simulation of roll response. The standard 
deviation of predicted roll response, 5.30, compared well with the measured roll 
response which is 5.29. A comparison of spectral densities, Figures 3.16, and 
histograms, Figures 3.17a and 3.17b, confirm that measured and predicted roll 
response compare quite well. Note that the there is a significant increase in 
probability mass at the extreme values. This behavior is typical of nonlinear systems 
with a softening stiffness (which is the case for barge roll motion). 

3.4 Form 2: Relative Motion Neglected in Damping Moment 

The second form of damping representation to be examined is the following 

! ui * * C 44JV*I*I *■ (2.S) 

- o 

In this form, the relative motion between the wave slope velocity and the barge roll 
velocity is neglected. As a result, only the barge roll velocity contributes to the 
equation. The following results will show this model displays an equal ability to 



roll (degrees) roll (degrees) 
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(a) Measured Wave Profile 
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(c) Measured Roll 




0 200 400 600 800 1000 

(d) Predicted Roll time (seconds) 



Figure 3.14 Comparison of Measured and Filtered Wave, Measured and Predicted 
Roll, White Noise Spectrum, Hs= 6 Ft (SB33, eqn 2.2) 
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Figure 3.15 Comparison of Spectral Densities (FC2/Hz), Measured and Filtered 
Wave, White Noise Spectrum, Hs= 6 Ft (SB33, eqn 2.2), - measured, — predicted 
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Figure 3.16 Comparison of Spectral Densities (Deg^2/Hz), Measured and Predicted 
Roll, White Noise Spectrum, Hs= 6 Ft (SB33, eqn 2.2), - measured, — predicted 



occurances occurances 
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(a) Measured Roll Amplitude (degrees) 



predicted roll 




(b) Predicted Roll Amplitude (degrees) 

Figure 3.17 Histograms of Measured and Predicted Roll, White Noise Spectrum 
(SB33, eqn 2.2) 



predict response amplitude and a marked improvement in matching the phase of the 
measured roll response. Table 3.2 lists the resulting damping parameters for this 
model. 
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Test Case 


Linear Damping 
Parameter 


Nonlinear Damping 
Parameter 


H=6.6 ft, T=5 sec 
(SB26) 


0.09 


0.07 


H=6 ft, T=6 sec 
(SB27) 


0.10 


0.07 


H=7.2 ft, T=8 sec 
(SB29) 


0.33 


0.40 


H=5.7 ft, T=10 sec 
(SB30) 


0.03 


0.07 


Bretschneider Spectrum 
(SB25) 


0.04 


0.04 


White Noise Spectrum 
(SB33) 


0.07 


0.08 



Table 3.2 Linear and Nonlinear Damping Parameters for Form 2, Equation 2.8 

3.4.1 Regular Wave, H= 6.6 Ft, T= 5 Sec (SB26) 

Figure 3.18 displays the results of the model barge subjected to a 6.6 ft wave 
with a 5 sec period considering only barge roll velocity. The damping parameters 
used to achieve these results are 0.09 for the linear term and 0.07 for the nonlinear 
term. The predicted amplitude closely matches the measured data. In addition, the 
predicted roll response motion is nearly identical in phase to that of the measured 



data. 
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measured and predicted wave and roll , H=6.6 ft , T=5 sec, SB26 




Figure 3.18 Comparison of Measured and Filtered Wave, Measured and Predicted 
Roll, Regular Wave, H= 6.6 Ft, T= 5 Sec (SB26, eqn 2.8) 

- measured wave, -. predicted wave, — measured roll, .. predicted roll. 
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3.4.2 Regular Wave, H= 6 Ft, T= 6 Sec (SB27) 

Figure 3.19 represents the barge response when subjected to a 6 ft, 6 sec wave 
considering only barge roll velocity in the damping term. The damping parameters 
used to achieve these results are 0.10 for the linear term and 0.07 for the nonlinear 
term. The response amplitude adequately matches the measured response of 11.5 
degrees. However, the predicted roll response is approximately 0.5 sec (1/12 period) 
out of phase with the measured data. Although not an exact match, this is an 
improvement over the phase difference of approximately 1 sec (1/6 period) produced 
in the relative motion damping form. 

3.4.3 Regular Wave, H= 7.2 Ft, H= 8 Sec (SB29) 

Figure 3.20 represents the barge response when subjected to a 7.2 ft, 8 sec 
period wave considering only the barge roll velocity in the damping term for 
predicted response. The damping parameters used to match the measured roll 
response of 4. 1 degrees are 0.33 for the linear coefficient and 0.40 for the nonlinear 
coefficient. As with the relative motion damping case, these coefficients are 
significantly higher than the previous test cases. Further study should be conducted to 
evaluate this phenomenon. This form of the model is an improvement over the 
relative motion damping form because it more adequately matches the measured 
behavior, particularly in the phase (at the peaks). 
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time (seconds) 

Figure 3. 19 Comparison of Measured and Filtered Wave, Measured and Predicted 
Roll, Regular Wave, H= 6 Ft, T= 6 Sec (SB27, eqn 2.8) 

- measured wave, predicted wave, -- measured roll, .. predicted roll. 
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Figure 3.20 Comparison of Measured and Filtered Wave, Measured and Predicted 
Roll, Regular Wave, H= 7.2 Ft, T= 8 Sec (SB29, eqn 2.8) 

- measured wave, -. predicted wave, — measured roll, .. predicted roll. 
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3.4.4 Regular Wave, H= 5.7 Ft, T= 10 Sec (SB30) 

Figure 3.21 represents the roll response of the barge subjected to 5.7 ft regular 
waves with a 10 sec period, considering only barge roll velocity in the damping term 
of the predicted response. As with the relative motion case, it is difficult to match 
response amplitude with the superharmonic effect. A damping parameter of 0.03 for 
the linear damping term and 0.07 for the nonlinear term provides the closest match to 
the measured roll response of 2 degrees. The standard deviation of the measured data 
is 1.84 compared with 1.87 for the simulated roll response. This is a small 
improvement over the prediction using relative motion damping which resulted in a 
standard deviation of 1.88. Figure 3.22 shows that the spectral densities of the 
measured and simulated roll response compare favorably for this form of the model. 

3.4.5 Bretschneider Spectrum, Hs= 6 Ft, Tp= 8 Sec (SB25) 

Figures 3.23a - 3.23d display the measured and simulated wave excitation and 
measured and predicted roll response. The damping parameters used in the model are 
0.04 for both the linear and nonlinear damping terms. Figures 3.23c and 3.23d 
indicate a good visual comparison between the measured and simulated responses. 

The standard deviation of 3.96 for measured data compared well to a predicted 
standard deviation of 3.96. The energy under measured and simulated spectral 
density curves are a close match although there is a slight shift in peak frequencies, as 
shown in Figure 3.24. A final comparison of probability densities, Figures 3.25a and 
3.25b, confirm that this form of the model provides a good simulation for the random 



wave environment. 
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Figure 3.21 Comparison of Measured and Filtered Wave, Measured and Predicted 
Roll, Regular Wave, H= 5.7 Ft, T= 10 Sec (SB30, eqn 2.8) 

- measured wave, -. predicted wave, - measured roll, .. predicted roll. 
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Figure 3.22 Comparison of Spectral Densities (Deg A 2/Hz), Measured and Predicted 
Roll, H= 5.7 Ft, T= 10 Sec (SB30, eqn 2.8), - measured, -- predicted 
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(a) Measured Wave Profile 




(b) Filtered Wave Profile 




0 200 400 600 800 1000 

(c) Measured Roll 




0 200 400 600 800 1000 



(d) Predicted Roll time (seconds) 

Figure 3.23 Comparison of Measured and Predicted Wave and Roll, Bretschneider 
Spectrum, Hs= 6 Ft, Tp= 8 Sec (SB25, eqn 2.8) 
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Figure 3.24 Comparison of Spectral Densities (Deg^2/Hz), Measured and Predicted 
Roll, Bretschneider Spectrum, Hs= 6 Ft, Tp= 8 Sec (SB25, eqn 2.8) 

- measured, — predicted 



occurances occurances 
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(a) Measured Roll Amplitude (degrees) 



predicted roll 




Figure 3.25 Histograms of Measured and Predicted Roll, Bretschneider Spectrum 
(SB25, eqn 2.8) 
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3.5.6 White Noise Excitation, Hs= 6 Ft (SB33) 

Figures 3.26a - 3.26d display the measured and filtered wave and predicted roll 
response for a white noise wave excitation with wave periods of 4 to 20 secs and a 
wave height of 6 ft. The damping parameters used in this model are 0.07 for the 
linear term and 0.08 for the nonlinear term. The resulting standard deviation for 
predicted roll response is 5.30 which is very close to the measured response standard 
deviation of 5.29. A comparison of spectral densities, Figure 3.27, and probability 
densities, Figures 3.28a and 3.28b, indicate the model provides an accurate prediction 
of roll response. Again a significant build-up in the extreme values is observed. 

3.5 Alternative Forms Considered 

Other forms of the damping moment studied are those represented by equations 
2.9 and 2.10: 

lji> + l m (4>~drildy) + C ul 4> + C MN (4>-dij/dy)\4>-drj/dy\ 9 ) 

+ R u (<l>,v,dv/dy) = 0 

+ -Wdy) + C UN (4> - di]fdy) \ 4> -dyldy | (2 .10) 

+ R u (<l>,V,Wty) = 0 

These forms are capable of simulating the response amplitude, however they are not 
as accurate in duplicating the measured phase. Figure 3.29a - 3.29d is a comparison 
of all four forms for test SB26. The two forms which use relative motion in the 
damping term have the greatest phase error. This indicates the system is inertia 
dominated and that linear (structural) damping makes a significant contribution to the 
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(b) Predicted Wave Profile 




0 200 400 600 800 1000 



(c) Measured Roll 




0 200 400 600 800 1000 



(d) Predicted Roll time (seconds) 



Figure 3.26 Comparison of Measured and Filtered Wave, Measured and Predicted 
Roll, White Noise Spectrum, Hs= 6 Ft (SB33, eqn 2.8) 
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Figure 3.27 Comparison of Spectral Densities (Deg A 2/Hz), Measured and Predicted 
Roll, White Noise Spectrum, Hs= 6 Ft (SB33, eqn 2.8)), - measured, — predicted 



occurances occurances 
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Figure 3.28 Histograms of Measured and Predicted Roll, White Noise Spectrum 
(SB33, eqn 2.2) 



roll (degrees) roll (degrees) r °M (degrees) roll (degrees) 
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(b) Form 2, Equation 2.8 




(c) Form 3, Equation 2.9 




(d) Form 4, Equation 2.10 time (seconds) 

Figure 3.29 Comparison of Roll Responses Using Four Different Forms (SB26) 

- measured, - predicted 
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response. The form which uses linear damping with roll velocity only and relative 
motion Morison nonlinear damping provides an improvement in the phase error, but 
is still slightly out of phase. Figure 3.29b represents the model which uses strictly 
structural damping. It is clear this model provides the best match. This trend is 
evident for the 6 and 8 sec wave cases as well as shown in Figures 3.30 and 3.31. 

The 10 sec wave test case displayed the same phase shifted results for all model forms 
as shown in Figure 3.32. 

In the random wave test cases, each model is equally successful at simulating 
the measured response. Equation 2.10, which uses strictly Morison damping, 
required a significant increase in the damping parameter for both random wave cases. 
For example with test SB25, the Bretschneider wave excitation, the damping 
parameters used for equation 2.2, 2.8, and 2.9 are very close to 0.04 for both the 
linear and nonlinear damping terms. The damping parameter for equation 2.10, 
which used only Morison damping, is 0.42. Similar results are obtained with test 
SB33, the white noise excitation. For equations 2.2, 2.8 and 2.9, the damping 
parameters are in the range of 0.07 to 0.08, where as the damping parameter for 
model equation 2.10 is 0.56. Tables 3.3 and 3.4 summarize the resulting damping 
parameters for equations 2.9 and 2.10. 

Form 4 of the SDOF model considered only relative motion Morison damping. 
It should be noted that all the damping parameters listed in Table 3.4 are extremely 
high for the regular wave cases. These ratios are unrealistic and indicate that Form 4 
is not an acceptable form of the model. 
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(b) Form 2, Equation 2.8 




(c) Form 3, Equation 2.9 




Figure 3.30 Comparison of Roll Responses Using Four Different Forms (SB27) 
- measured, — predicted 
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(a) Form 1, Equation 2.2 





(c) Form 3, Equation 2.9 




(d) Form 4, Equation 2.10 time (seconds) 

Figure 3.31 Comparison of Roll Responses Using Four Different Forms (SB29) 

- measured, -- predicted 
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(a) Form 1, Equation 2.2 




(b) Form 2, Equation 2.8 
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(c) Form 3, Equation 2.9 




(d) Form 4, Equation 2.10 time (seconds) 

Figure 3.32 Comparison of Roll Responses Using Four Different Forms (SB30) 

- measured, - predicted 
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Wave Test 


Linear Damping 
Parameter 


Nonlinear Damping 
Parameter 


H=6.6 ft, T= 5 sec 
(SB26) 


0.09 


0.14 


H=6 ft, T=6 sec 
(SB27) 


0.10 


0.11 


H=7.2 ft, T=8 sec 
(SB29) 


0.34 


0.34 


H=5.7 ft, T=10 sec 
(SB30) 


0.03 


0.07 


Bretschneider Spectrum 
(SB25) 


0.04 


0.04 


White Noise Spectrum 
(SB33) 


0.07 


0.08 


Table 3.3 Linear and Nonlinear Damping Parameters for Form 3, Equation 2.9 



Wave Test 


Nonlinear Damping 
Parameter 


H=6.6 ft, T=5 sec 
(SB26) 


1.07 


H=6 ft, T=6 sec 
(SB27) 


1.25 


H=7.2 ft, T=8 sec 
(SB29) 


88.00 


H=5.7 ft, T=10 sec 
(SB30) 


1.25 


Bretschneider Spectrum 
(SB25) 


0.42 


White Noise Spectrum 
(SB33) 


0.56 



Table 3.4 Linear and Nonlinear Damping Parameters for Form 4, Equation 2.10 
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4.0 Sensitivity Studies 

4.1 Variation of Damping Parameters of Form 2, Equation 2.8 

A sensitivity study was performed on the model which provided the best results 
for roll response, equation 2.8, in which damping from the wave excitation was 
disregarded and only structural damping was considered. For each wave test, the 
nonlinear damping parameter is held fixed and the linear damping parameter varied by 
a wide range to determine the effects on simulated response. The linear parameter is 
then held fixed and the nonlinear damping parameter is varied by a wide range to 
determine the resulting effects on simulated roll response. 



4.1.1 Regular Wave, H= 6 Ft, T= 5 Sec (SB26) 

The previously established coefficients for this case were a linear damping 
parameter of 0.09 and a nonlinear damping parameter of 0.07. Table 4.1. shows a 
constant decrease in standard deviation as damping parameters increase for both the 
linear and nonlinear variations. 



Linear 

Damping 

Parameter 


Nonlinear 

Damping 

Parameter 


Standard 

Deviation 


Linear 

Damping 

Parameter 


Nonlinear 

Damping 

Parameter 


Standard 

Deviation 


0.03 


0.07 


9.89 


0.09 


0.01 


9.54 


0.05 


0.07 


9.77 


0.09 


0.03 


9.48 


0.07 


0.07 


9.59 


0.09 


0.05 


9.42 


0.09 


0.07 


9.33 


0.09 


0.07 


9.33 


0.11 


0.07 


9.00 


0.09 


0.09 


9.28 


0.13 


0.07 


8.60 


0.09 


0.11 


9.20 


| 0.15 


0.07 


8.15 


0.09 


0.13 


9.12 



Table 4.1 Variation in Damping Parameters and Resulting Standard Deviations (SB26) 
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Figures 4. la - 4. le display the amplitude and phase changes which result from 
varying the linear damping parameter. The trend observed from Figure 4. 1 is that 
with the lower linear damping parameter, the response amplitude is greater than the 
measured. In addition, the phase of the predicted response lags the measured by 
approximately 0.25 secs, Figure 4.1a. As the linear damping parameter proceeds 
through the range, the amplitude decreases and the phase shift transitions from a 0.25 
sec lag to leading the measured response by approximately 0.25 secs, Figure 4. le. 

The effects of the variation of the nonlinear damping parameters is not as 
evident. Figures 4.2a - 4.2e display the change in roll response with variation in the 
nonlinear damping term. The predicted response phase is almost identical to the 
measured response. There is not a visual change in the amplitude as the nonlinear 
parameters proceed through the range. However, a comparison of standard deviations 
shows a slight decrease as the nonlinear parameter increases, indicating this term does 
provide a contribution to the response. 

4.1.2 Regular Wave, H= 6Ft,T=6 Sec (SB27) 

The linear and nonlinear damping parameters for this case were established as 
0.10 and 0.07 respectively. Table 4.2 indicates a constant decrease in standard 
deviation as the damping parameters increase. Figure 4.3 represents the resulting 
changes when the nonlinear damping parameters were held fixed and the linear 
damping parameters were varied. It is evident that amplitude decreases with an 
increase in damping parameter. Figure 4.4a - 4.4e represents the results of holding 
the linear damping parameter fixed and varying the nonlinear damping parameter. 

The simulated response amplitude decreased through the range of damping 



roll (degrees) roll (degrees) roll (degrees) roll (degrees) roll (degrees) 
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parameters. It appears the phase shift remains unchanged throughout the range. 



Linear 

Damping 

Parameter 


Nonlinear 

Damping 

Parameter 


Standard 

Deviation 


Linear 

Damping 

Parameter 


Nonlinear 

Damping 

Parameter 


Standard N 
Deviation i 


0.04 


0.07 


9.82 


0.10 


0.01 


8.76 


0.06 


0.07 


9.41 


0.10 


0.03 


8.67 


0.08 


0.07 


8.97 


0.10 


0.05 


8.60 


0.10 


0.07 


8.50 


0.10 


0.07 


8.50 


0.12 


0.07 


7.90 


0.10 


0.09 


8.43 


0.14 


0.07 


7.59 


0.10 


0.»1 


8.33 


0.16 

1 


0.07 


7.47 


0.10 


0.13 


8.25 



Table 4.2 Variation in Damping Parameters and Resulting Standard Deviations (SB27) 
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4.1.3 Regular Wave, H= 7.2 Ft, T= 8 Sec (SB29) 

The damping parameters established for test case SB29 are 0.33 and 0.40. 
There is a constant decrease in standard deviation as the damping parameters increase 
as shown in Table 4.3. 



Linear 

Damping 

Parameter 


Nonlinear 

Damping 

Parameter 


Standard 

Deviation 


Linear 

Damping 

Parameter 


Nonlinear 

Damping 

Parameter 


Standard 

Deviation 


0.20 


0.40 


3.40 


0.33 


0.25 


3.06 


0.25 


0.40 


3.26 


0.33 


0.30 


3.05 


0.30 


0.40 


3.12 


0.33 


0.3t 


3.04 


0.33 


0.40 


3.03 


0.33 


0.40 


3.03 


0.40 


0.40 


2.84 


0.33 


0.45 


3.03 


0.45 


0.40 


2.70 


0.33 


0.50 


3.02 


0.50 


0.40 


2.56 


| 0.33 


0.55 


3.01 



Table 4.3 Variation in Damping Parameters and Resulting Standard Deviations (SB29) 



Predicted response amplitude is greater than measured when the nonlinear 
damping parameter is held fixed and the linear parameter is varied. The phase of the 
predicted response leads the measured at the lower linear damping parameters. As 
the linear parameter proceeds through the range, the phase shift changes and 
eventually the predicted response lags the measured by approximately 0.5 sec. These 
results are displayed in Figures 4.5a - 4. 5e. Holding the linear damping parameters 
fixed and varying the nonlinear produces no significant change in amplitude or phase 
difference, as shown in Figures 4.6a - 4. 6e. 
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4.1.4 Regular Wave, H= 5.7 Ft, T= 10 Sec (SB30) 

Test case SB30 is the only test that displays sensitivity to change in damping 
parameters. Table 4.4 shows the change in standard deviation over the range of 
damping parameters. The sensitivity is very evident in Figures 4.7a - 4.7e which 



Linear 

Damping 

Parameter 


Nonlinear 

Damping 

Parameter 


Standard 

Deviation 


Linear 

Damping 

Parameter 


Nonlinear 

Damping 

Parameter 


Standard 

Deviation 


0.01 


0.07 


2.55 


0.03 


0.01 


1.89 


0.02 


0.07 


2.054 


0.03 


0.03 


1.88 


0.03 


0.07 


1.87 


0.03 


0.05 


1.87 


0.05 


0.07 


1.74 


0.03 


0.07 


1.87 


0.07 


0.07 


1.70 


0.03 


0.09 


1.86 


0.09 


0.07 


1.68 


0.03 


0.11 


1.85 


0.11 


0.07 


1.66 


0.03 


0.13 


1.84 



Table 4.4 Variation in Damping Parameters and Resulting Standard Deviations (SB30) 



show the rapid change in simulated response. The superharmonic virtually disappears 
once the linear damping parameter reaches 0.05. Figure 4.8 displays the entire 1024 
sec time series for this case. It is evident that the damping parameter has a major 
effect on the occurrence and magnitude of the superharmonic. The standard deviation 
changes rapidly with the increased damping parameter until the superharmonic is 
eliminated at 0.05 damping. This is emphasized by the fact that when the linear 
damping parameter varies from 0.01 to 0.05, the difference in standard deviation is 
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Figure 4.7 Variation of Linear Damping Parameter, Constant Nonlinear Damping 
Parameter, £ N = 0.07, H= 5.7 Ft, T= 10 Sec (SB30, eqn 2.8) 

— measured, .. predicted 
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0.8025. However, when the damping parameter is varied from 0.05 to 0.11 the 
change in standard deviation is only 0.08. A variation in the nonlinear damping 
parameter has very little effect on the response amplitude or superharmonic response 
as shown in Figure 4.9. 



4.1.5 Random Wave Tests 

Tables 4.5 and 4.6 show the sensitivity of the model to Bretschneider and white 
noise wave excitation. Both damping models exhibit a significant change in standard 
deviation over the range of linear damping parameters. In each case, the standard 
deviation decreases by over 50% as the linear parameter increases. 



Linear 

Damping 

Parameter 


Nonlinear 

Damping 

Parameter 


Standard 

Deviation 


Linear 

Damping 

Parameter 


Nonlinear 

Damping 

Parameter 


Standard 

Deviation 


0.01 


0.04 


6.15 


0.04 


0.01 


4.14 


0.02 


0.04 


5.16 


0.04 


0.02 


4.07 


0.03 


0.04 


4.42 


0.04 


0.03 


4.01 


0.04 


0.04 


3.96 


0.04 


0.04 


3.96 


0.06 


0.04 


3.35 


0.04 


0.06 


3.86 


0.08 


0.04 


2.97 


0.04 


0.08 


3.77 


0.10 


0.04 


2.70 


0.04 


0.10 


3.69 



Table 4.5 Variation in Damping Parameters and Resulting Standard Deviations (SB25) 



71 




a) | N = 0.01 




b) In = 0.03 




c) In = 0-07 





e) In = 0.13 time (seconds) 



Figure 4.9 Variation of Nonlinear Damping Parameter, Constant Linear Damping 
Parameter, | L = 0.03, H= 5.7 Ft, T= 10 Sec (SB30, eqn 2.8) 

- measured, .. predicted 
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Linear 
H Damping 
Parameter 


Nonlinear 

Damping 

Parameter 


Standard 1 
Deviation 


Linear 

Damping 

Parameter 


Nonlinear 

Damping 

Parameter 


Standard 

Deviation 


0.01 


0.08 


8.45 


0.07 


0.02 


5.67 


0.03 


0.08 


6.98 


0.07 


0.04 


5.53 


0.05 


0.08 


6.01 


0.07 


0.06 


5.41 


0.07 


0.08 


5.30 


0.07 


0.08 


5.30 


0.09 


0.08 


4.74 


0.07 


0.10 


5.20 


0.11 


0.08 


4.29 


0.07 


0.12 


5.10 


0.13 


0.08 


3.93 


0.07 


0.14 


5.01 



Table 4.6 Variation in Damping Parameters and Resulting Standard Deviations (SB33) 



4.2 Variation of Nonlinear Restoring Moment 

All of the models considered in this analysis have used a 13th-order polynomial 
restoring moment. A sensitivity test is conducted to determine the effects on the 
predicted roll response when lower-order polynomials are used to fit the actual 
restoring moment curve. 



4.2.1 Cubic Fit to Restoring Moment Curve 

Figure 4. 10 shows the actual restoring moment curve and the 3rd-order 
polynomial used to provide an analytical fit to the curve. This cubic term is then 
substituted into equation 2.8. The resulting predicted roll response is shown in Figure 
4.11. A comparison of the measured and predicted spectral densities for roll 
response. Figure 4.12, indicates this polynomial is acceptable although it does not 
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x 10 s 




Roll (deg) 



Figure 4.10 Comparison of Numerical (Paulling, 1995) and Analytical (Yim et al. 
1995) Restoring Moment Curves, Cubic Fit, - numerical, -- analytical. B 13 5 = 
3.9502e6, B n>5 = -5.1930e6 
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Figure 4.11 Comparison of Measured and Predicted Roll, Cubic Restoring Moment 
H= 7.2 Ft, T= 8 Sec (SB29), - measured, .. predicted 




frequency (Hz) 

Figure 4.12 Comparison of Spectral Densities (Deg^2/Hz), Measured and Predicted 
Roll, Cubic Restoring Moment, H= 7.2 Ft, T= 8 Sec (SB29) 

- measured, — predicted 
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provide and "exact" match of the magnitude of the predict roll response relative to the 
ones with higher-order terms(see following section). 

4.2.2 5th, 9th and 13th-Order Restoring Moment Curves 

Figure 4.13 shows the 5th-order polynomial fit to the restoring moment curve. 
The resulting predicted roll response is shown in Figure 4. 14. A comparison of the 
measured and simulated spectral densities, Figure 4.15, confirms that the 5th-order 
restoring moment polynomial provides accurate results. Figures 4. 16 and 4.17 show 
the 9th and 13th-order polynomial fits to the restoring moment curve and the resulting 
roll responses, shown in Figures 4.18 and 4.19. 

This sensitivity study indicates that a 3rd or 5th-order polynomial fits the 
restoring moment curve with sufficient accuracy to provide valid roll response 
simulation. However, it should be noted that it does not require a significant amount 
of time or effort to use the 13th-order fit, and this high-order polynomial becomes a 
requirement to accurately predict roll response at large roll angles. 



Righting Moment (ft-lb) 
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x 10 s 




Roll (deg) 



Figure 4.13 Comparison of Numerical (Paulling, 1995) and Analytical (Yim et al. 
1995) Restoring Moment Curves, 5th Order Restoring Moment, 

- numerical, — analytical. B 135 = 5.5384e6, B n5 = -1.1747e7, B 95 = 5.2136e6 
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Figure 4.14 Comparison of Measured and Predicted Roll, 5th Order Restoring 
Moment, H= 7.2 Ft, T= 8 Secs (SB29), - measured roll, .. predicted roll 




Figure 4.15 Comparison of Spectral Densities (Deg”*2/Hz) , Measured and Predicted 
Roll, 5th Order Restoring Moment, H= 7.2 Ft, T= 8 Sec (SB29) 

- measured, — predicted 
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Figure 4.16 9th Order Fit to Restoring Moment Curve. - measured, — actual 
B 13>5 = 5.7927e6, B n>5 = -1.4128e6, B 9>5 = 1.0693e6, B^ -3.9235e6, 

B 5i5 = 5.5127e6 




Ron (deg) 

Figure 4.17 13th Order Fit to Restoring Moment Curve. - numerical, — analytical 
B U5 = 5.6051e6, B u5 = -1.16577, B<, >5 = 1.7258e6, Bj S = 9.6915e6, 

B 5 , 5 = -9.3298e6, B 3>5 = 3.3999e6, B 15 = -4.4569e5 
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Figure 4. 18 Comparison of Measured and Predicted Roll, 9th Order Restoring 
Moment, H= 7.2 Ft, T= 8 Secs (SB 29), -- measured roll, .. predicted roll 




Figure 4.19 Comparison of Measured and Predicted Roll using 13th Order Restoring 
Moment, H= 7.2 Ft, T= 8 Secs (SB29), -- measured roll, .. predicted roll 
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5.0 Comparison of Single, 2- and 3- Degree-of-Freedom Models 

The final objective of this study is to compare the results of roll response for a 
3DOF and a 2DOF model with the SDOF model. A comparison of the 5, 8 and 10 
sec regular wave test cases is done as well as the Bretschneider spectrum. Table 5.1 
lists a comparison of damping parameters for each of these cases. It is evident from 



Wave 

Test 


SDOF 

Linear 

Damp- 

ing 

Parame * 
ter 


SDOF 

Non 

linear 

Damp- 

ing 

Parame 

ter 


2DOF 

Linear 

Damp- 

ing 

Parame 


2DOF 

Non 

linear 

Damp- 

ing 


3DOF 

Linear 

Damp- 

ing 

Parame 


3DOF 
Non 
linear 
Damp- 
ing 1 


ter 


Parame 
ter 1 


ter 


6 

Parame 

ter 


SB26 


0.09 


0.07 


0.32 


0.32 1 


0.40 


0.40 


SB29 


0.33 


0.40 


0.55 


0.55 


0.25 


0.25 


I SB30 


0.03 


0.07 


0.02 


0.01 


0.03 


0.03 


| SB25 


0.04 


0.04 


0.08 


0.08 


0.05 


0.05 



Table 5.1 Damping Parameters for SDOF, 2DOF and 3DOF Models 



this table that there is a wide range of damping parameters among the different 
models. There does not appear to be any trend established by the damping 
parameters for each wave test. In some instances, the parameters are higher in the 
SDOF model and in others the 2DOF and 3DOF exhibit higher parameters. 

Figures 5.1a - 5.1c, which show response to 5 sec wave excitation, indicate 
that all three models simulate the amplitude of roll response very well. The 3DOF 
and SDOF model match the phase identically, however, the 2DOF model leads the 
measured response slightly. 
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Figure 5.1 Comparison of 3DOF, 2DOF and SDOF Models, H= 6.6 Ft, T= 5 Sec 
(SB26) 
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All three models accurately predict the roll response for the 8 sec wave 
excitation case, Figures 5.2a - 5.2c. The amplitude as well as phase are nearly 
identical in each case. It is interesting to note that each of the models required 
unusually high damping parameters for this particular wave case. 

Figures 5.3a and 5.3b compare SDOF and 3DOF responses for the 10 sec wave 
test case. This is the only case observed where the higher-degree-of-freedom model 
exhibits better prediction capabilities than the SDOF one. It should be noted that the 
3DOF and 2DOF models used analytical input and the SDOF used measured input. 
The SDOF model predicts the subharmonic activity seen in Figure 5.4b. Figure 5.4 
compares the entire time series of the 2DOF and SDOF models. Again, the multi- 
degree-of-freedom model predicts superharmonic response much more accurately. 
However, results are not available to show prediction capabilities using measured 
input for all three models. 

The final comparison is for the Bretschneider Spectrum test case. This is the 
only test case in which all models exhibit damping parameters which are close in 
magnitude. Figures 5.5a - 5.5c indicate the SDOF and 3DOF models provide an 
accurate prediction in this wave case. The 2DOF model, Figure 5.5b, is observed to 
be less capable of producing the same results. 

A more thorough study is recommended to determine the adequacy of each 
model and to perhaps determine some correlations between the models. These 
preliminary results, however, indicate the SDOF model may be equally as accurate 
and in some cases more accurate than the multi-degrees-of-freedom models. Use of 
the single degree of freedom model can be an excellent tool for future research. 




(a) 3D0F Model, predicted results 
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(c) SDOF Model time (seconds) 

Figure 5.2 Comparison of 3DOF, 2DOF and SDOF Models, H= 7.2 Ft, T= 8 Sec 
(SB29) 
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(a) 3DOF Model, predicted results 




(b) SDOF Model 



Figure 5.3 Comparison of 3DOF and SDOF Models, H= 5.7 Ft, T= 10 Sec 
(SB30) 
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(b) SDOF Model 



Figure 5.4 Comparison of 2DOF and SDOF Freedom Models, H= 5.7 Ft, T= 10 
Sec (SB30) 
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(a) 3DOF Model, predicted results 





(c) SDOF Model time (seconds) 

Figure 5.5 Comparison of 3DOF, 2DOF and SDOF Models, Bretschneider Spectrum 
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6.0 Concluding Remarks 

In this report, four different forms, each identical except for a variation in the 
damping term were studied. From these equations, an SDOF model to predict barge 
roll motion was identified. The form which provided the most accurate predicted 
responses proved very capable of simulating barge response in both the regular and 
random wave environment with the exception of test SB30. A sensitivity study was 
conducted to determined the effects of varying damping parameters. In addition, the 
restoring moment term was examined to determine the effects of using 3rd, 5th, 9th, 
and 13th-order terms. Finally, a comparison was made between the SDOF, 2DOF 
and 3DOF models to identify the strengths and weaknesses of each model among the 
different wave test cases. 

6.1 Conclusions 

1) The model which provided the most accurate response simulation was 
equation 2.8. This form assumes the relative motion effects can be neglected. As a 
result, only the barge roll velocity is considered in the damping terms. 

2) The SDOF model was capable of predicting barge response for all regular 
and random wave test cases except SB30, in which the 10 sec wave excitation was a 
multiple of the barge’s natural roll frequency. As a result, the response motion 
contained a superharmonic which was difficult to simulate. 

3) A constant damping parameter cannot be used to predict all responses 
throughout the wave excitation spectrum. The damping parameters for the regular 
and random wave excitations summarized in Table 6.1, show no obvious pattern of 



behavior. 
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Wave Test 


Linear Damping 
Parameter 


Nonlinear Damping 
Parameter 


H= 6.6 Ft, T= 5 Sec 
(SB26) 


0.09 


0.07 


H= 6 Ft, T= 6 Sec 
(SB27) 


0.10 


0.07 


H= 7.2 Ft, T= 8 Sec 
(SB29) 


0.33 


0.40 


H= 5.7 Ft, T= 10 Sec 
(SB30) 


0.03 


0.07 


Bretschneider Spectrum 


0.04 


0.04 


White Noise Spectrum 


0.07 


0.08 



Table 6.1 Summary of Damping Parameters for Form 2, Equation 2.8 



4) A sensitivity study of the restoring moment term determined that a 3rd or 
5th-order polynomial is sufficient for simulation when the magnitude of the roll 
response is not large, although the 3rd-order polynomial tends to slightly underpredict 
response at times. There is no discemable difference in predicted response provided 
by a 9th or 13th-order restoring moment polynomial. The incremental computational 
time using these higher-order polynomials is minimal. In addition, when roll 
response becomes large, the higher-order polynomials become a requirement. 

5) When comparing the SDOF model with the multi-degree-of-freedom 
models, the results are mixed. In the regular wave tests, the SDOF model provides 
the same results as the 3DOF, except in the 10 second wave test case. The ability of 
SDOF and 3DOF models to simulate random response is remarkably better than the 
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2D0F model. These observations combined with the significant time saving makes 
the SDOF very desirable to use. 

6.2 Recommendations for Future Studies 

1) Further studies in variation of the damping terms may improve the ability 
to predict superharmonic response. Varying the effects of relative motion, roll 
velocity and wave slope velocity could possibly produce better results. 

2) A computer program is currently under development to predict the 
probability of barge capsize. A comparison of the SDOF, 2DOF and 3DOF models 
can be used in this program to establish further strengths and weaknesses of each 
model. Although simplistic compared to the 2DOF and 3DOF models, the SDOF 
model may provide equally adequate results with much less time invested. 

3) In determining the damping parameters in this study, the standard deviation 
of the amplitude was used as the criterion for determination. This assumes the phase 
shift is unimportant. The focus of an alternate study could be to use a least-square 
criterion which weighs both amplitude and phase shift in the error calculation. 

4) In the random wave cases, a comparison of spectral densities showed a 
shift in peak frequencies between the measured and predicted spectrum. Variations in 
input data such as dominant period of wave spectrum could be analyzed to determine 
if this shift may be eliminated or decreased. 
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